Fluctuational transitions through a fractal basin boundary 
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Fluctuational transitions between two co-existing chaotic attractors, separated by a fractal basin 
boundary, are studied in a discrete dynamical system. It is shown that the mechanism for such 
transitions is determined by a hierarchy of homoclinic points. The most probable escape path from 
the chaotic attractor to the fractal boundary is found using both statistical analyses of fluctuational 
trajectories and the Hamiltonian theory of fluctuations. 
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The mechanism of fluctuational escape from a 
chaotic attractor (CA) through a fractal basin 
boundary (FBB) represents one of the most chal- 
lenging; unsolved problems in fluctuation theory 
HI Q Q 01 The unpredictable and highly com- 
plex stochastic behavior of such systems arises in 
part from the presence of limit sets of complex 
geometrical structure, and in part from the frac- 
tality of the basin boundary [fj El ■ For this rea- 
son, the central question - whether or not there 
exists a generic mechanism for fluctuational tran- 
sitions through the FBB - has remained unan- 
swered. More specifically, it has remained unclear: 
(i) if boundary conditions can be found both on the 
CA and on the FBB; (ii) if there exits a unique es- 
cape path from the CA to the FBB; (iii) whether 
this path can be determined using the Hamiltonian 
theory of fluctuations; (iv) if there is any determin- 
istic structure involved in the transition through 
the FBB itself; and (v) what effect is exerted by the 
noise intensity. If transitions across FBBs are char- 
acterised by general features, a knowledge of them 
could considerably simplify analyses of both sta- 
bility and control for chaotic dynamical systems, 
which areproblems of broad interdisciplinary in- 
terest Hit 

A promising approach to the solution of this 
problem is based on the analysis of fluctuations 
in the limit of very small noise intensity. In this 
limit, a stochastic dynamical system fluctuates to 
remote states along certain most probable deter- 
ministic paths 
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, corresponding to rays in 
the WKB-like asymptotic solution of the Fokker- 
Planck equation [l2j. The possibility of extending 
such an approach to chaotic systems, both continu- 
ous and discrete, was established earlier 
It was shown also that the presence of homoclinic 
tangencies, causing the fractalization of the basins, 
causes a decrease in the activation energy [l3j. 

In this Letter we show that a generic mechanism 
of fluctuational transition between co-existing CAs 



separated by an FBB does exist, that it is deter- 
mined by a hierarchy of homoclinic original saddles 
forming the homoclinic structure, and that there is 
a unique most probable escape path (MPEP) from 
the CA that approaches an accessible orbit on the 
fractal boundary. 

To demonstrate the existence of this escape 
mechanism, we take as an example the two- 
dimensional map introduced by Holmes The 
properties of this map, including the structures 
both of its CA and of its locally disconnected FBB, 
are generic for a wide class of maps and flow sys- 
tems EH Ell- It is this fact, taken with the results 
of our investigations of escape in other systems, 
that allows us to conclude that the escape mecha- 
nism described below is indeed a typical one. The 
Holmes map is 

Xn+l = Vn (1) 

Vn+i = -bx n + dy n -yl + £ n , 

where is zero-mean, white, Gaussian noise of 
variance D. Due to symmetry, the noise- free sys- 
tem (JTJ with b = 0.2 and 2.0 < d < 2.745 has 
pairs of co-existing attractors, the basins of which 
are separated by a boundary that may be either 
smooth or fractal, depending on the choice of pa- 
rameter values. We choose for our studies b = 0.2 
and d — 2.7, which corresponds to there being 
two co-existing CAs separated by a locally discon- 
nected FBB (see Fig. [TJ. The fractal dimension 
of the boundary has been determined numerically 
(dim = 1.84472) by using the "uncertainty expo- 
nent" technique introduced in E3- The chaotic 
attractors in JQ) appear as the result of a period- 
doubling cascade and, for the parameters chosen, 
each of them consists of two disconnected parts. 

We have modelled the system Q numerically, 
exciting it with weak noise, and have collected tra- 
jectories that include escape paths from one CA 
to the other, and also the corresponding realisa- 
tions of noise that induced these transitions. By 
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FIG. 1: The co-existing chaotic attractors (filled black 
regions) and their basins of attraction represented by 
grey and white respectively. The accessible boundary 
saddle points of period 3 are shown by the small filled 
black circles labelled S3. Their stable manifolds are 
shown as full black lines. The saddle points of period 
1 are shown by the crosses labelled SI. The saddle 
point at the origin is labelled O. 

ensemble averaging a few hundred such escape tra- 
jectories and noise realisations, we have obtained 
the optimal escape path and the corresponding op- 
timal force, which are shown in Fig. [5] The results 
of this statistical analysis allow us both to deter- 
mine the boundary conditions near the CA and 
the FBB, and to demonstrate the uniqueness of the 
MPEP. It can be seen in particular that, in leaving 
the CA, the system Q falls into a small neighbour- 
hood of the saddle point of period 1 (SI) located 
between its two disconnected parts and having the 
multipliers p 1 = 0.118975 and p 2 = 1.681025. Its 
stable manifolds separate the parts of the CA, 
while the unstable ones belong to the CA. The 
system makes a few iterations in the neighbour- 
hood of SI (initial plateau in Fig.|5fa)) and then 
moves to the FBB in three steps, crossing it at 
a saddle point of period 3 (S3) with multipliers 
pi = 0.001016 and p 2 = 7.875512. Calculations 
have shown that, for the chosen parameter values, 
S3 lies within the FBB. Moreover, its stable mani- 
fold (solid black line in Fig.QJ is dense in the FBB 
and detaches the open neighborhood, including an 
attractor, from the FBB itself. This allows us to 
classify it as an accessible boundary point 0] . 

An analysis of the structure of escape paths in- 
side the FBB has shown that the homoclinic saddle 
points play a key role in its formation. In the sys- 
tem (QJ, we observe an infinite sequence of saddle- 
node bifurcations of period 3, 4, 5, 6, 7..., which oc- 



cur at parameter values dz < d\ < d$ < d§ < dj... 
and are caused by tangencies of the stable and un- 
stable manifolds of the saddle point O at the ori- 
gin. The homoclinic orbits appearing as a result 
of these bifurcations were classified earlier as orig- 
inal saddles, and it was also shown that their sta- 
ble and unstable manifolds cross each other in a 
hierarchical sequence ^| . To characterize this hi- 
erarchical relation between original saddles, it is 
reasonable to introduce a parameter p equal to 
the ratio | ln(pi(P)) | /\n(p 2 (P)), where p\ and 
p 2 are the multipliers of a saddle point P. Cal- 
culations have shown that, for the original sad- 
dles of period 3,4,5,6,7,8... in QJ, the follow- 
ing hierarchical sequence of index p values oc- 
curs: pz — 3.339, P4 — 3.08, p$ — 2.999, pa = 
2.339, fi 7 = 1-958, fig = 1.539. Moreover, the val- 
ues of index p corresponding to the other homo- 
clinic saddle cycles are close to zero. Correspond- 
ingly the probability of finding the system in their 
neighbourhood tends to zero. 

These results allows us to infer features of a fluc- 
tuational transition through a locally disconnected 
FBB that are probably generic, as follows: (i) it al- 
ways occurs through a unique accessible boundary 
point; and (ii) the original saddles forming the ho- 
moclinic structure of the system play a key role 
in the formation of the paths inside the FBB, the 
difference in their local stability defining the hier- 
archical relationship between them. Thus, we may 
claim that complicated structure of escape trajec- 
tories, caused by the thin homoclinic structure and 
their randomness, has in many respects a deter- 
ministic nature. 

Having now understood the mechanism of es- 
cape, we can seek the MPEP. According to the 
Hamiltonian theory of fluctuations Q, 0, 0, the 
MPEP is the path which minimizes the energy 

n=l 

of the possible realizations of noise inducing 
a transition of the system from the CA (with 
the initial condition on SI) to the FBB (with the 
final condition on the accessible orbit S3). The La- 
grangian of the corresponding variational problem 
can be found following (cf. 0) in the form 

N N 
n—1 n— 1 

where Q is taken into account using the Lagrange 
multiplier A„. Varying L with respect to A„, 



and x n , the following area-preserving map is ob- 
tained: 

Vn+i = -bx n + dy n -yl + X y n (3) 
= (d-3xl +1 )K/b~X»Jb 

A n+1 — \i 

Equations are supplemented by the following 
boundary conditions 

lim A£ = 0, (x n ,y° n ) e SI, (x^yl) € S3. (4) 

n — > — oo 

The MPEP is the minimum-energy heteroclinic 
trajectory connecting SI to S3 in the phase space 
of @ . The solution of this boundary value prob- 
lem is in general complicated, because of the pres- 
ence of multiple local energy minima [2(ij | induced 
by the complex geometrical structure of the unsta- 
ble manifolds of S3 (see e.g. jij for a discussion). 
The solution of the boundary value problem in- 
volves a parameterization of the structure of the 
multiple local minima requiring, in turn, a proper 
parameterization of the unstable manifold in the 
vicinity of the initial conditions 22\. The MPEP 
found by this method is shown in Fig. [21 It can 
be seen that, within the range shown by the verti- 
cal dotted lines in Fig.l^a), the theoretical MPEP 
closely coincides with the path obtained by statis- 
tical analysis of escape trajectories in the Monte 
Carlo simulations. Note that no further action is 
required to bring the system to the other attrac- 
tor once it has reached the accessible orbit of the 
FBB, i.e. once it has reached the points numbered 
8 in Fig. 12 a) and (b); correspondingly, the opti- 
mal force measured in the numerical simulations 
(inset) falls back to zero. 

The existence of an almost deterministic mecha- 
nism of transition across the FBB raises important 
questions about the effect of noise on this mecha- 
nism, and on the structure of escape trajectories 
inside the FBB. We have therefore used randomly 
chosen initial conditions in a very small neighbor- 
hood of the accessible point S3 through which es- 
cape occurs (see Fig. 2(b)). By definition, any ar- 
bitrarily small neighborhood of S3, lies within the 
FBB, and must contain points belonging to the 
basins of both attractors. Therefore the system 
can cross the FBB starting from a very small neigh- 
borhood of S3, even in the absence of noise. By 
collecting all such successful escape paths, we have 
calculated the probabilities for the system to pass 
via small neighborhoods of different original sad- 
dle cycles during its escape, both in the presence 
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FIG. 2: (a) the most probable escape path (dashed 
line) connecting the CA with the period-3 saddle cy- 
cle lying on the fractal boundary, obtained from the 
Monte-Carlo simulations with D = 10~ J . The opti- 
mal path found by the solution of the boundary-value 
problem is shown as a solid line. The ^-coordinate of 
the saddle point SI is shown by the horizontal dashed 
line, (b) A two-dimensional plot of the paths pre- 
sented in (a) where the results obtained by solution of 
the the boundary value problem (consecutively num- 
bered points indicated by circles, corresponding to the 
numbered points in (a)) coincide almost perfectly with 
those obtained by numerical simulation (stars). Inset 
in (a): the optimal force as determined in the numeri- 
cal simulation. 



and absence of noise. As can be seen from Fig. |3| 
the corresponding probabilities demonstrate the 
same hierarchical interrelationship in both cases, 
which is determined by the value of index /j, de- 
fined above. This structure is evidently robust 
with respect to noise-induced perturbations. The 
addition of noise causes a slight broadening of the 
distribution in Fig. [3] and a small increase of the 
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FIG. 3: Probabilities of finding a fragment correspond- 
ing to the different period-T original saddle cycle in the 
collected escape trajectories. 



probability for the system to escape via original 
saddles of larger period. 

In conclusion, we have described the mechanism 
by which noise-induced escape occurs through a 
locally disconnected FBB. We have found the 
(unique) most probable escape path from a chaotic 
attractor to the fractal boundary, using both sta- 
tistical analyses of fluctuational trajectories, and 
the Hamiltonian theory of fluctuations. We have 
shown that the original saddles forming the homo- 
clinic structure play a key role in effecting the tran- 
sition through the FBB itself. In particular, their 
local stability defines the hierarchical relationship 
between the probabilities for the system to pass 
via small neighborhoods of different original sad- 
dle cycles during its escape, both in the presence 
and in the absence of noise. We emphasize that the 
escape mechanism we have revealed must be appli- 
cable to the broad class of two dimensional maps 
and flows 0, 0, 0] that exhibit the same type of 
FBB. For instance, one possible application of our 
results is to the development of an energy-optimal 
control scheme for the CO2 laser, a discrete model 
of which demonstrates the type of FBB considered 
above @. 
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